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Possibility to model the high energy synchrotron emission (in X- and 7-rays) 
' from supernova remnants (SNRs) is an important task for modern astronomy and 
astrophysics, because it may be responsible for the nonthermal X-rays and TeV 7- 

■ rays observed recently from a number of SNRs. This emission allows as to look in 
the processes of particle acceleration on SNR shocks and generation of cosmic rays. 
In this paper, a model for the synchrotron emission from shell SNR in nonuniform 
interstellar medium and nonuniform magnetic field is presented. This model is a 
generalization of the model of Reynolds & Chevalier developed for a spherical SNR 
in the uniform medium and uniform magnetic field. The model will be used for 
studies on the thermal and nonthermal X-ray images and spectra from nonspherical 
SNRs in different interstellar magnetic field configurations. 
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^ 1 Introduction 

For the last few years, a considerable interest has grown in the investigation of the 
synchrotron X-ray emission from supernova remnants (SNRs), caused by electrons 
accelerated on SNR shocks up to the energy ~ 10 13-14 eV. This emission process 
is believed to be responsible for the nonthermal X-ray (SN 1006, Cas A, Tycho, 
G347.3-0.5, IC443, G266.2-1.2, RCW86) and TeV 7-ray (SN 1006, Cas A, G347.3- 
0.5) emission observed recently from a number of SNRs. 

Model for synchrotron X-rays from shell SNRs developed up to now [TTL H works 
with uniform magnetic field and is based on the simple spherical Sedov model 
for SNR hydrodynamics. The most essential simplifications of the Sedov model are 
the uniform interstellar medium and spherical SN explosion. 

Observations certify that SNRs are nonspherical objects with complicated dis- 
tributions of surface brightness. The nonuniform interstellar medium are the most 
important factor in modification of the structure and evolution of SNRs. We present 
here a model for the synchrotron emission from shell SNR in nonuniform medium 
and nonuniform magnetic field. This model is a generalization of the model Jl|, [2| . 

2 Model for synchrotron emission from nonspher- 
ical SNR in nonuniform magnetic field 
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2.1 Hydrodynamics 



The position R(9,p) and velocity D(9,ip) of the shock front as well as the distri- 
bution of gas parameters inside the volume of SNR depend on the ambient density 
distribution p°(r,9,(p) and anisotropy of the supernova explosion energy E sn (9,(p). 
It is enough to have parameters of gas within R(9, p) calculated to model the ther- 
mal X-ray emission of SNR. Besides p°(r, 9, p) and E sn (9, (p), the synchrotron image 
and spectrum of SNR are also affected by the distribution of interstellar magnetic 
field B D (r, 9, ip) which may also be nonuniform. 

It is assumed hereafter that SNR are on the adiabatic phase of evolution. The 
dynamics of the shock front of nonspherical adiabatic SNR and the distribution of gas 
parameters inside whole shocked volume can be described with the hydrodynamic 
method presented in 0. In this method, a three-dimensional (3-D) object is divided 
on a number of 1-D sectors and the distributions of parameters are found separately 
for each of them. Note, that the velocity of flow is radial in the hydrodynamic 
method M and the direction of a sector and electron velocity in it are given by the 
same angles (9, <p) of a spherical coordinate system with the origin in the place of 
explosion. 

We assume also that 7 = 5/3, shock compression ratio a = 4, gyrofactor 77 = 1 , 
Bqmb/B — > 0, maximum energy attained before the adiabatic stage -E m ax,o = (see 
for details). We do not consider the effect of the cosmic ray pressure what can 
change a . 

2.2 Emission 

The synchrotron volume emissivity of a fluid element (r, 9, p) is || 

S v = C{a)KB^ +1),2 v-^l\ (1) 

where C(a) is a constant, B± is the tangential component of the magnetic field 
(perpendicular to the electron velocity in the fluid element), v is a frequency, K is 
the normalization of the sharply truncated electron distribution 

f KE~ a dE, E < E max 
N(E)dE = „ max . (2) 

I U, ±j > -C/ max 

Let us consider the distribution of the emission parameters over the shell and 
their evolution downstream. 



2.2.1 Spectral index 

The power index a is constant downstream because electrons lose energy proportion- 
ally to their energy (Eq. fli~0|) ) and remain essentially confined to the fluid element 
in which they were produced. The power index a is also constant over the surface 
of a nonspherical SNR because the shock is strong. Namely, in the first order Fermi 
acceleration mechanism, a = (a + 2)/(er — 1), where the shock compression ratio o 
does not depend on the ambient density distribution in the strong shock limit (Mach 
number 3> 1) realized in SNRs: o = p s /p° — (7 + l)/(7 — 1); the subscripts "s" and 
"o" refer to the values at the shock and to the surrounding medium respectively. 
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2.2.2 Magnetic field components 



The components of the magnetic field evolve differently behind the shock: B± rises 
everywhere at the shock by the same factor o = p s / p° and B\\ is not modified by the 
shock g : 

B ± , a = B° ±>s -^ o , B U = B° U . (3) 

r s 

No further turbulent amplification of the magnetic field is assumed. 

Downstream variation of the components follows from the flux-freezing condition 
B^(r)rdr = B^_(a)ada, continuity equation p(r)r 2 dr = p°(a)a 2 da and magnetic flux 
conservation r 2 B\\(r) = a 2 B°(a), where a is Lagrangian and r = r(a,t) is Eulerian 
coordinates. Thus, if both B and p° are nonuniform, the components of magnetic 
field B are 

BxM = B±M eM ^ BfM = w (4) 

The components of the ambient magnetic field B Q in the point (r, 6, if) are 

B° ± = B \sine \, Bl = S |cos6 |, (5) 

where the obliquity angle B Q is the angle between the directions of the magnetic 
field B D in this point and the electron velocity in the sector. 

2.2.3 Normalization K 

The assumption that the energy density to of relativistic particles is proportional to 
the energy density of the magnetic field in each fluid element 



E 

UJ = 

R ■ R 



J max ^max 

J EN(E)dE = K J E 1 ~ s dE oc B 2 (6) 



gives for the variation of K over the shell and downstream 



K oc < (7) 
B 2 \n(E min /E max ), a = 2. 



2.2.4 Injection energy 

Particles start to accelerate on the shock when their energy E become more than 
the injection energy E m i n . The first order Fermi acceleration mechanism suggests for 
E m in- E min = 1.68 10~ 4 In A 2 / 3 (X n tot /D) 2 ^ eV, where In A is Coulomb logarithm, 
A G = 1.0 10 18 Tg /n tot cm is the mean free path of electron, n tot is the total number 
density (electrons plus ions) p. This E min is 14.65 In A 2 / 3 times the average kinetic 
energy of electrons. 

Variation of E m \ n over the shell of nonspherical SNR finds for each sector inde- 
pendently: 

E min (R, 6, ip) = 0.026 In A 2/3 D 3 (6, f) 2 MeV, (8) 
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where D 3 is the shock velocity in 10 3 km/s. 

If particles move downstream leaving the region of acceleration, the energy varies 
as§ 

E(a) =E -p(a) l /\ (9) 
if electron loses its energy due to the adiabatic expansion only 

E = E 7>/(3p) (10) 

where p(r) = p(r)/p s , E D is the energy of the fluid element a at time t Q when the 
shock crosses it. 



2.2.5 Maximum energy 

The maximum energy E max which electrons are accelerated to varies over the shell 
of nonspherical SNR. Therefore -E max also finds independently in each sector, with 
the expressions given in f2]| : E max (R, 9, <p) — min(.£/ maXi i, E max ^2, -Emax, 3 

) where the 

radiative losses, finite time of acceleration and escaping of electrons allow for the 
maximum energies: 

E maxA = 7.810 13 F 1 B- 1 J 2 D 3 eV, (11) 
£ max , 2 = 1-9 10 14 f F 2 B ^ Dj dt 3 eV, (12) 

-'tad 

E maX) 3 = 7.7 10 12 A maXj i7 B 0tlM eV. (13) 
Here the age of SNR t and the time of transition into the adiabatic stage t a d(#, <p) 

'4' 



are in 10 3 years, B Q41 is B Q in pG, A maXj i 7 is A max in 10 17 cm. This A max ~ 10 16 18 cm 



is a free parameter. 

The geometric factors are || 



with 



F^JG/Rj, F 3 = l/Rj, (14) 



Z + A/a B = cos 2 9 + 1 
Z + Aa B ' cos 2 9 + r 

_ l ( T i? (cOS 2 e o +l)+4(cOS 2 + l) 

Rj = 2 4T^ ' (15) 



B s l + 16tan 2 o 



B° V l + tan 2 6 ' 

where the obliquity angles are 9 D for upstream and 6 for downstream; tan 
4tan0 . 

The maximum energy changes downstream in accordance to (|9|) . 
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2.2.6 Obliquity angle 



Synchrotron emission of the fluid element (r, 9, (p) depend on the values B a (r, 9, 4>) 
and o (r, #,</>). The calculation of O depend on the way how the distribution 
B D (r, 9, cf>) is given. 

If we know distributions of B Q (r, 9, (ft), 9 H (r, 9, (p) and <fH( r , 8, 0), where (9s, ipn) 
are the spherical angles of the vector B D in the point (r, 9, 0), then 

cos O = sin 9 sin 9 H cos(y? — ipu) + | cos 9 cos 9 H \ . (16) 

There are three independent of a coordinate system angles in our task: the 
inclination angle 8 between the line of sight and the density gradient, the aspect 
angle between the line of sight and the ambient magnetic field, and the third 
angle £ between the magnetic field and the density gradient. 

Let as consider two interesting particular cases. 

i) . If direction of the ambient magnetic field is uniform, i.e. 0(r, 0, 0) = const, 
then let us accept that the origin of the coordinates coincides with the place of 
explosion, the axis x of Cartesian coordinate system is oriented opposite to the 
direction of the line of sight (toward observer) and that B Q lies in the plane (xz). 
In such a case cos O is given by (|1^) with 9h = vr/2 — and cpu = 0. 

ii) . If direction of the density gradient is uniform, i.e. 8(r, 9, 0) = const, then let 
us accept that the origin of the coordinates coincides with the place of explosion, 
axis z of Cartesian coordinate system is oriented opposite to the density gradient 
and the line of sight lies in the plane (xz). Geometrical consideration yields for 
8 ^ and £ ^ 

cos0 o = sin 9 sin £ cos(ip — ifu) + | cos 9 cos £|, 

. cos0 — I cos 5 cos £| (17) 

where cos if h = ; r • 

sin 8 sin £ 

If 8 = or 7r (gradp° is parallel to the line of sight) but £ ^ 0, then 6 H = n — £ 
and — £ respectively, and additionaly ^(r, 9, (p) should be known. If £ = or tt 
(direction of B Q coincides with gradp°), then O = 7r — 9 or 9 respectively. 



3 Conclusions 

The model for the synchrotron emission of nonspherical SNR in nonuniform inter- 
stellar magnetic field is presented. Being applied to studies of high energy emission 
of SNRs, in particular to X-rays, it allows as to analyse the thermal and synchrotron 
X-ray images and spectra of the objects, to make conclusions about the SNR itself, 
supernova explosion, structure and properties of the interstellar medium, shock wave 
physics, acceleration processes on the shocks, 7-ray and cosmic ray generation. 
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